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1. Introduction 

In the present day the cosmology faces two important problems. The first is that the 
common matter cannot explain the galaxy rotation curves and so another type of matter 
is needed pQ. Although this is an old problem [2], the nature of this another type of 
matter is unknown until the present and several candidate particles are proposed in 
the literature [3]. Nowadays the most accepted explanation that cosmologists have for 
the galaxy rotation curves is that there exists a strange matter field that interacts only 
gravitationally with the known matter - the so-called dark matter. Some alternative 
theories to the dark matter fluid are the Modified Newtonian Dynamics (MOND) [3], 
the f(R) [5] and fi(R) [6] (non-minimally coupled matter) theories which consider the 
dark matter as a geometrical effect. The recent data from the measures of the matter 
contained in the galaxies through the gravitational lensing effect strongly support the 
existence of the dark matter [7J [8] . 

The second problem is that the astronomical observations show that the expansion 
of the Universe is accelerated in the late time [91 [10] . The standard cosmology cannot 
explain what causes this acceleration and cosmologists have to look for a satisfactory 
explanation for this observed fact. The first attempt of explaining it came from 
the revival of the cosmo logical constant. Unfortunately, such an idea soon presented 
inconsistences [H] and new models appeared in the literature, but no one has proved 
to be the definitive one yet [12]. Among these models, a well known is the quintessence 
model [131 E] > which can present several acceptable forms for the potentials, proposed 
phenomenologically or inspired by fundamental theories [151 [HI ELY]- The models with 
fermion fields [HI 13 [20], tachyon fields [211 122] (k-essence), Chaplygin gas [231 [23] 
and van der Waals gas [251 EE] are also alternatives to the cosmological constant. All 
these models essentially try to model an exotic component with negative pressure which 
composes the most part of the energy density of the Universe - the so-called dark 
energy. The f(R) theories [271 [28] also propose an explanation of geometric origin for 
the accelerated expansion. 

We starts this work from an action with boson (quintessence-type) and Dirac 
fermion fields minimally coupled to the gravity and a common matter field. One 
considers a spatially flat, homogeneous and isotropic Universe (in according to the 
observations) and a point-like Lagrangian is obtained from the action. The undefined 
potentials in the Lagrangian are selected from the existence condition for the Noether 
symmetry - interesting results using this approach can be found in the references 
[221 EH [301 EU [321 EH [34]. The selected potentials indicate that the boson field can 
describe the dark energy and an additional matter field and the fermion field behaves 
as the standard matter. From this result we describe the dark sector by bosons and 
fermions, with the bosons composing a dark energy-matter field and the fermions 
playing the role of dark matter. This type of fermionic dark matter (but with additional 
interactions) is also studied in [35]. The standard matter behavior of Dirac fermion fields 
can be seen in other contexts [321 ESI [37j . 
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Since the potentials satisfying the Noether condition imply that the Lagrangian 
presents a Noether symmetry, there is a constant of motion and a cyclic variable can 
be obtained through a coordinate transformation in the configuration space. This can 
help us integrate the dynamical system. Then, using these additional informations 
furnished by the Noether symmetry, the field equations of our original problem are 
completely integrated. The result is that the Universe described by the proposed 
action - satisfying such a symmetry - evolves with alternated periods of accelerated 
and decelerated expansion. 

The organization of the paper is the following: in the second section the action 
of the model is presented and the field equations from a point-like Lagrangian for 
a flat Friedmann-Robertson- Walker (FRW) metric are derived. The selection of the 
potentials from the Noether symmetry condition is performed in the third section and 
the identification of the fermion field with the dark matter is discussed in the fourth 
section. In the fifth section the field equations are integrated via constant of motion 
and cyclic variable and the analysis of the corresponding solution is fulfilled in the sixth 
section. The paper is closed with the conclusions in the seventh section. 

In this work we adopt the signature (+, — , — , — ) and the natural units 87rG = c = 
h = l. 



2. Action and field equations 

Let us start from a general action describing a boson field and a fermion field minimally 
coupled to the gravity 



S 




:<T^<R0 - U(4>) 



i 

+ 2 



M 



(1) 



where represents a boson field, ip and ip = ^7° denote the spinor field and its adjoint, 
respectively, which represent a Dirac fermion field, and Cm describes a common matter 
field. R is the Ricci scalar, the dagger denotes complex conjugation transpose and 
$ = ipif). Further, U and V are the self-interaction potentials of the boson and fermion 
fields, respectively. In (JTJ the covariant derivatives and _D M are defined by 

= (<9 M + V, D$ = d$ + #V (2) 

where is the spin connection 

1 



(3) 



with denoting the "vierbein" and T M = e^7 a . The 7° are the Pauli-Dirac matrices. 

By varying the action with respect to the metric tensor g^ u one obtains the 
Einstein's field equations 



Rfj,v 2 9nvR 



-T, 



flUl 



(4) 
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with = T*5 + T$ v + T* v , where the letters M, ip and label the energy- momentum 
tensors of the common matter, boson and fermion fields, respectively, which are defined 
by 



1 V» t /=H hn»v ' v>> 



T£, = 7 [ #V£>^ + #VA,V> - - A^VV 1 - &hA> ( 6 ) 
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T*, = ^<R0 - - C/) ^, (7) 

where £^ represents the Dirac Lagrangian 

c f = % - [ ^d m v - (PrfW^ ] - v. (8) 

From the action (JTJ, for a flat FRW metric, with the fields <fi and ip spatially 
homogeneous and a pressureless common matter field, we can write a point-like 
Lagrangian of the form 

£ = 3aa 2 - a 3 - U^j + % -o? (^ T > - ^ 7 °V>) + a 3 V + p° M , (9) 

which furnishes the same dynamical equations that those from the Einstein's field 
equations The point represents time derivatives and a is the scale factor. Here 
p° M is the energy density of the common matter at an initial instant (this comes from 
the common matter energy density: pu = Pm/ ?)- 

Imposing that the "energy function" associated with the Lagrangian (jHJ) is null, the 
result is the Friedmann's equation, i.e., 

dC. dC- ^dC dC ; . n 

E c = jrra + -jcj) + ^— + — - C = 0, 

oa d(f) dtp dip 

a\ 2 _ Pm + Pj, + P<t> ^ 



with and p^ denoting the energy densities, defined as 

Pi> = V 1 P4 ,= l -tf + U. (11) 

In according to the energy-momentum tensors and ([7]). 

From the Euler-Lagrange equation for a, applied to fl9]), one has the acceleration 
equation 

a _ Pm + Pj> + P<$> + 3 (P± + P<t>) ( 12 ) 
a 6 

and for ip and ip, we have the Dirac's equations for the spinor field and its adjoint 
coupled to the gravitational field, respectively, 

* + l(!)* + ^ = °- < 13 > 

— 3 /d\ — — n dV . . 
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In (I12p the pressures are denned as 

Also in according to (jBJ) and (j7j). 

Finally, the Euler-Lagrange equation for gives 

/d\ • dU , , 

^ +3 O + ^ =0 - (16) 

which is just the Klein-Gordon equation for the field 0. 
3. Selection of potentials 

One can express the Lagrangian in terms of the components of the spinor field 

ip = (tpx, -02, -03, ip4) T and its adjoint = (0>*,0>2; ~i J h ~ t PD- So it takes the form 



C = 3ad 2 - a 3 ( %- - U J + V £ (0*0, - 0*^) + a¥ + p° M . (17) 



Thus the dynamics of the system is now described by 10 coordinates, i.e., the 
configuration space of the system is represented by {a, 0, 0>* , ipA, with j = 1, 2, 3, 4. 

The undefined forms of the potentials can be selected by the requirement of the 
existence of a Noether symmetry for a given Lagrangian. This approach constrains 
the possible functions for the potentials at the same time that it provides a conserved 
quantity associated with the selected forms. Since the dynamical system presents a 
Noether symmetry, one has a constant of motion and a cyclic variable is allowed, which 
can help us integrate the system. A Noether symmetry exists for a given Lagrangian 
of the form C = £(qk,<jk) if the condition L^C = is satisfied, with L x designating 
the Lie derivative with respect to the vector field X (called infinitesimal generator of 
symmetry) defined by 

d dau d \ , . 



X 



dq k dt dq k ) 

where the a^s are functions of the generalized coordinates q k . 
From the E-L equations one arrives at the relation 

!(? a '3 =i * £ = o < (19) 

and the conserved quantity (or constant of motion) associated with the Noether 
symmetry is obtained 

dC 

M » = X>%- (20) 
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The new set of variables {Qk(Qi)} f° r the configuration space, such that one of the 
variables is cyclic, obeys the following system of differential equations 



(21) 



where z is the cyclic variable. The range of the index I is the same of that of k and the 
Uk-i and z comprehend the new coordinates Qk( a i) of the configuration space. 

Starting from the Noether symmetry condition L^_C = X£ = applied to the 
Lagrangian (ITT)) , with the infinitesimal generator of symmetry X defined for our problem 

as 



8 • d 8 ■ d 

8a 8a oq) d(j) 

4 / 8 8 -8 ■ 8 \ 

fr[\ dip? 8^ dipt dipij 

we have the following system of partial differential equations 

C + 2a "^- = °> 3C + 2a ~Q0- = °> 6 ^ ~ a = °> I 23 ) 



A/<9C <9A,*\ n A/ 



4 /r)C- r)D \ 

3Co'0j + aD, + a g (^ - ^,*J = 0, (25) 

3CWf + ( ^- ( ,g(^-^)=0, (26) 

gc aco aA ^A f27 x 

3C (U + V) + aD ^ + a£ + D t e^*) ^ = 0, (28) 



where the symbol e« takes the values 



6j = +1 for 2 = 1,2; 

ei = -l for z = 3,4. (29) 

One can see from equations ()2Tj) that the coefficients Co and A are functions only 
of a and <fi. Then, assuming that Cq and A are separable functions 



C = ci(a)c 2 (0), A = d 1 (a)d 2 ((f}), 



(30) 
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we obtain the solution for the system 



Co = , (31) 



D <> = - J ■ (32) 
Cj = -a ^ + 06^*, (33) 

A e a<f> _|_ £> e -act> 



[/ = u Q (Ae^ - Be' 01 *) , V = V V, (35) 



with A, B, Uq, Vq and /3 being constants and a — 



2 V 2 



These are the coefficients and potentials that satisfy the Noether symmetry 
condition. Now, if the potentials (|35|) are taken for the Lagrangian (|T7|) . it presents 
a Noether symmetry. 



4. Dark matter as a fermion field 



The solution for the potential of the fermion field is V = Vq^, which is essentially a term 
of mass. Thus we can replace the constant Vo for m in order to identify this constant 
with the mass of the fermions. Once it is done, we determine from equations ( TTTT) i and 
f fT5|) i that the fermion field has its energy density and pressure given by 

Pi) = m ty, pip = 0, (36) 

characterizing a field of pressureless matter. Note that this form for requires that ^ 
must be a positive quantity. 

Observe that we did not assume a priori that the fermion field is a massive field, 
whereas the potential takes this information. This characteristic of the field comes 
exclusively from the symmetry condition when it constrains the potential to a term of 
mass. 

This result suggests that the fermion field behaves as a standard matter field. But 
obviously this field has a nature different from that of the common matter since it is 
composed exclusively by fermionic particles. Thus, having in view that these fermions 
produce an additional pressureless matter, one may identify the fermion field with the 
dark matter. For this propose, we must assume that this field describes fermionic 
particles that interact only gravitationally with the common matter or have very weak 
non-gravitational interactions with it. 

On the other hand, the field <fi received from the Noether condition the potential 
( l35]) i . which allows an accelerated expansion. Hence, from this point on, the dark sector 
will be identified with the fields and ip. 
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5. Searching for exact solutions 



Taking U and V given by ( 1351) . we have to solve the following system of coupled 
differential equations 

2 2 

U Q (Ae a * - Be-**) , (37) 



»( 


d\ 


2 






a/ 




a 3 








2 










2° 


+ 


(3 


i + 


a 









^2 



2 u 

Y -^{Ae^ - Be' a(p y = 0, (38) 

+ 3 + 2af/ (A 2 e 2a ^ - 5V 2q *) = 0, (39) 

ip + - (-} ip + irnfip = 0, % + - (-} ~ip - irmprf = 0. (40) 
2 \aj 2 \a/ 

Once with these forms for U and V the dynamical system presents a Noether 
symmetry, one has an additional dynamical equation provided by the constant of motion, 
namely 

dC 



_ dC <9£ ( „ dC 
M = C — + -D —j + £ Ci-f- + A- 



<9Vv 

= (Ae a * + £e _a *) v/oad + (Ae a(t> - Be -a *) a 3/2 + ip a 3 *, (41) 

which is determined from equation (1201) . In this equation one can make /3 = z/3o, with /3o 
being a real constant, in order to have a constant of motion whose value is real, whereas 
\l> is a real positive quantity. 

As it was exposed before, if the system presents a Noether symmetry, there is 
a cyclic variable that reduces the dynamics of the system, which is obtained from a 
transformation of variables. The transformation that takes the configuration space from 
{a, <j>,ipj ,ipj} to {z,u,Vj,Wj} satisfies the system 

^P + D ^ + Hc^ + Di ^-)=0, (42) 
da d<P £i\ d^* d^J 

^ + A^ + Ef^ + A^)=0, (43) 
da d(j) dipt d^J 

+ DM + ± ( • D*f) 0, (44) 
da d(j) dip* diPiJ 



da 9^ til dipt dtp J 



(45) 



with z being the cyclic variable. This system is obtained by applying f[2~Tj) . 

Using the constant of motion (l4T|) and a transformation of variables satisfying (142!) - 
15]), we will look for a solution to the field equations (!37|) -( l40|) . 
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Firstly, the Dirac's equations (jlD|) 1 and flUUh can reduced to an equation for \1/ 



* = 0, 



(46) 



which is easily integrated in terms of a, giving 

vj/ = — - 

a 6 

where ^ is a positive constant. Or, by solving equation (j40]) i in terms of a, the solution 
for ip is 



(47) 



( ^}\e- imt \ 

^0 e -tmt 



ip°e imt 



-3/2 



V fie imt j 



with being constants that satisfy the relation 



( 4 * 



In this point it is interesting to note that the solution \1/ = ^o/a? implies that the 
fermion field has energy density = mty /a 3 . So the extra matter field is naturally 
added with the common matter field. This first result corroborates to our identification 
of the fermion field with the dark matter. 

Having in view that the Dirac's equations could be solved independently of the 
others field equations, in terms of a, one has effectively a reduction of the system ( 1371) - 
( 140]) to differential equations involving only the dynamical variables a and <fi. Then we 
will solve the problem through the reduced system 



d\ 
aj 



p° M + mtf c 



+ 



a<f> 



+ 3 



^ + U (Ae 
<p + 2aU (A 2 e 2a * - B 2 e~ 2o " t ' 



Be 



-Ot(j> 



[Ae^ + Be~ a ^) VQa~d + (Ae^ - Be: 



0. 

3/2, 



M , 



where M = M + Ad^o and the result \l/ = ^o/a 3 was used. 

Hence it is enough to find a new set of variables related only to a and 
adequate particular solution of this type for the system ( T42|) -( l45j) is 

,Ae a<t> - Be'"* 



u 



,3/2; 



,3/2 



y/6AB 



0. 



W i 



0. 



(49) 

(50) 
(51) 

An 

(52) 
(53) 



V6AB 

Using ( |53|) we can express the constant of motion in terms of the cyclic variable 
obtaining 

. M 



4AB' 



(54) 
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This equation is immediately integrated, furnishing 

z(t) = z 1 t + z 2 , (55) 

with z\ = Mq/AAB and z 2 = constant. 

The Friedmann's equation takes the following form in the new variables 

•2 -2 , o A T3TT 2 , PO + m ^0 

z = u + 3ABU u H — — — . (56) 

2AB 

Substituting z from equation (1541) into the Friedmann's equation, we get its integration 
and obtain the solution 



u(t) = u sm (cut + b ), (57) 



where 



,2 _ M - SAB (p°, + mfp) n 2_o ARU 

u 48A z B 3 Uq > <" " ^^o, (58) 

and b is a constant. 

For the solutions z(t) and expressed in terms of the original variables through 
the relations fl52l) i and fl53|) . we easily get the explicit forms of a(t) and 0(t) 

f to 2 \ ^ 3 1/3 
a(t) = f — J {zjt 2 + 2z x z 2 t + z 2 2 -ul sin 2 (cot + 6 ) } , (59) 

0(t) = -1 In { ^ + ^ + u sin (ut + b ) 1 _ 1 ^ ,A y 



Finally, from equation ffjHj) one can write the time evolution of the spinor field 
\/2T7o 



( ip°e- tmt \ 

^O e -tmt 



i>le imt 



3 C 

where 

9(t) = {zjt 2 + 2 Zl z 2 t + z 2 2 -ul sin 2 (ut + 6 )} • (62) 
6. Constraining the solution 

Now let us analyze the behavior of a. The solution (|59j) is valid from the time when 
the regime is matter-dominated. Close to this time we have that the time evolution of 
a is near the form a cx t 2//3 , i.e., when the oscillatory term UqSui 2 (cot + bo) presents its 
minimum value and the term z\t dominates z 2 in the expression fl59|) . We can satisfy this 
limit situation by choosing the constant of integration z 2 = and replacing the constant 
bo for -ur, such that sin 2 [u(t — r)] has its minimum value (zero) at t = r, where r is 
the time when one has a matter-dominated regime. Before performing the analysis, one 



G(t), (61) 
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introduces the dimensionless variables: t = H t, f = H t, m = jjr, A = B = -J^, 



3ABU , u = H u , M 



Mo 



z 2 = H z 2 and Q,° M 



p° M , Wf 

where if is the Hubble parameter at the present time. Hence the solutions (J59"j) . (I5T71) 
and (THTj) take the forms 



a(t ) : 

0(t) 



UJ 



1/3 



2C/q, 



m sm \u 



(t-r)}} 



1/3 



1 j f zit + «o sin [co(t — f)] 
2a | z\t — uq sin [u(t — r)] 



2a 1 5. 



v 7 ^ 



3 

„4ini 



{^P-Uq sin 2 [w(t-f)]V 



-1/2 



(63) 
(64) 

(65) 



/ 



V tie 1 ' 

From (163]) . the normalization a(trj) = 1 implies the relation 
2 2f/ + MqW 2 sin [o)(t — r)] 



(66) 



with to being the present time. By substituting z\ = M /AAB into fl66l) and using fl58|) i 
in the dimensionless form, written as 







3<7 (3U M 2 - 24u 2 n° Mi) 
16w 6 



one obtains M in terms of A, B and Z7o, namely, 

2 8u 2 {Au 2 - 9Sl° m sin 2 [u(t- f )] } 



Mr 







9t/ {^o-sin 2 [o;(t-r)]} 



(67) 



(68) 



Thus we have three free parameters A, B and Uq that determine all the other 
parameters u, M , u and z\ since one knows f^^to and f. Using the solution (|M|) to 
rewrite the energy density and pressure of the field 0, we have the following dimensionless 
expressions for the energy densities and pressure 



Pu{t ) 



P°M 

H$a(t )3 
ui 4 z 2 ul 



Aa 2 Uia(t f 

2w 4 ugsin 2 [u(t- f)] 



ait) 3 ' rvvw H 2 a(t) 3 
{sin [u(t — t)\ — ut cos [u>(t— f)] } 



3^, 
a(t)3 : 

2 



(69) 



+ 



4 ) = P<t>(t ) - 



3U a(t ) 3 

4cD 4 ^ sin 2 [u(t- r)] 



(70) 



(71) 



3C/ a(t ) 3 

If we observe the solution (|63|) and the set of equations (169]) - (ITT]) and take into 
account (16"6"]) - (16"8"]) . we conclude that the time evolution of the quantities (scale factor, 




Figure 1. Density parameters of the matter fields (common matter and fermion field) 
and the boson field. Up to t = 1.5 (left frame) and t = 0.002 (right frame). 



energy densities and pressure) that will furnish the cosmological parameters depend 
essentially on Uq and the product AB. Once the Hubble time is defined as to — l/-Ho an d 
our dimensionless time is t — Hot, the approximated value for the dimensionless present 
time is to — 1. Note that the dimensionless present time is not exactly equal to the unity 
due to the fact that the Universe presents a non-constant rate of expansion (decelerated 
expansion when t < 1 and accelerated expansion when t > 1). By taking the Hubble 
time as about 14 Gyr (see e.g. p2J[38]) and the beginning of the matter-dominated era 
when the Universe is about 75 kyr old, f takes the value f = 5.4 x 10~ 6 to = 5.4 x 10~ 6 . 
The oscillating term of a has its period determined by u>, which can be constrained by 
adopting Q° M = 0.04, 0° = 0.22 and = 0.74. These values are based on the references 
[j~2| 139] such that at t we observe Q = 0.04 + 0.22 + 0.74 = 1 in our analysis. So, from 
these cosmological parameters valued at the present time, one can infer the values for 
the free parameters. In order to present the results for a realistic description - such that 
the observational data are respected -, the best values that we have found are Uq = 15 
and AB = 10~ 3 . 

In Fig. 1 are plotted the density parameters of the matter fields (comprehending 
the common matter and the fermion field) and the boson field as functions of t - which 
is just denoted by t from this point on. All possible values for Uq and AB imply that 
Q° M , and are reached at the times larger than t = 1. This shows that the solution 
naturally produces an accelerated Universe independently on the values of Uq and AB. 
Heaving in view that the age of the Universe is close to the Hubble time, the density 
parameters in Fig 1 are plotted in such a form that Q° M = 0.04, = 0.22 and = 0.74 
appear as close to t — 1 as it is possible. In Fig. 1 the present values of the density 
parameters are at t = 1.08 (~ 15 Gyr), which is 8 percent larger than the Hubble time. 
So in the following graphics the present time will be t = 1.08, characterizing a Universe 
with accelerated expansion. In the left of the Fig. 1 one can observe that is dominant 
today, but going to the past the matter fields become the dominant components. In the 
right of this figure we note that close to r one has that dominates the matter density 
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Figure 2. Density parameters of the boson field (left frame) and matter fields (right 
frame) for t > 1.5. 



parameters. However, from the solution of a we have that closing to r the approximate 
behavior is a oc t 2 ^ 3 , i.e., the Universe is matter-dominated. Then the boson field must 
play the role of a matter field at this time. This will be clarified by the next results. 

The density parameters up to the distant future of the Universe are plotted in Fig. 
2. In the left of this figure the evolution of is presented, which oscillates tending 
to a constant value for t — > oo. In the right frame the same type of evolution can be 
observed for QmiP, obviously. Let us observe that in the future fl^ will always dominate 

The time evolution of the equation of state of the boson field = p^/ptf, is plotted 
in Fig. 3. From the left frame we infer that evolves from a positive value and 
decreases to —1 and turns to increase slowly by reaching the value co^ = —0.73 at the 
present time - it is larger than that of the cosmological constant equation of state, 
uja = —1- Note that, according to [12], the value of the parameter inferred from the 
observational data for the dark energy has still great uncertainties. For times close to 
the beginning of the matter dominated era, one observes from this figure that going back 
to t the parameter reverses its sign and tends to 1. Thus it presents a behavior of 
matter, explaining the approximated form a oc t 2 ^ 3 for the scale factor with a dominant 
Vt^ in Fig. 1. In the right of this figure the time evolution of is presented for large 
times. From it one can infer that oscillates with time in the range — 1 < u,p < 1, 
which shows that the boson field behaves as matter and dark energy alternately. 

In Fig. 4 the acceleration of the scale factor is plotted. We started from the 
requirement that Q%,fl^ and must be observed near t — 1, being found precisely 
at t = 1.08, as it was allowed by the solution. Now, in the left of Fig 4, one can see 
that we naturally have an accelerated expansion that have passed through a transition 
from a decelerated to an accelerated expansion. This transition occurs about t = 0.51 
(7.1 Gyr) and the present period is accelerated. Going back to r, an decelerated period 
corresponding to a real matter dominated Universe occurs. The acceleration of the scale 
factor for the distant future is presented in the right frame. This curve shows a Universe 
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t t 

Figure 3. Equation of state of the boson field. Up to t = 1.5 and t — 0.002 (left 
frame), and for t > 1.5 (right frame). 




t t 

Figure 4. Acceleration of the scale factor. Up to t = 1.5 (left frame) and t = 50 (right 
frame). 

that will return to a decelerated period and turn to accelerate in the future and so on - 
this behavior is obviously caused by the oscillating equation of state of the boson field. 
The final result is that the solution describes a Universe that evolves with alternated 
periods of accelerated and decelerated expansion and this oscillatory behavior is a kind 
of deadened oscillation: when t — > oo, one has a — > 0. 

Now we can ask if going back in time the big bang nucleosynthesis is permissible, 
that is, if a radiation-dominated era can be obtained from this model. The solutions 
analyzed above showed that the field <fi plays the role of the dominant matter field in 
the matter-dominated era, with the equation of state of the form p^ = p^, informing 
us that at this time the referred field is kinetic dominant. To verify if a radiation 
era is possible from these results, we must solve the field equations for the boson and 
radiation fields. From what we know, such a system does not present an analytical 
solution for the potential U used here, then we will search for a numerical solution of this 
scenario. Firstly, we transform the derivatives with respect to time in the Friedmann, 
acceleration and Klein-Gordon equations into derivatives with respect to red-shift by 
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Figure 5. Density parameters of the boson (dashed line) and the radiation (straight 
line) fields as functions of z. 



using the relationships 

z 



1 

a 



7t = - H(l + z) Tz 



(72) 



and divide all of them by Hq. From this, through a little algebraic manipulation, we 
obtain the following system of coupled differential equations to solve 

2ffl'(l + z) = p R + Pi + p R + Pd>, (73) 



= PR + P<t>+PR + P</>, 

H 2 (l + z) 2 <P" + H{l + z) \H\l + z) 



2H 



,/ dU 

4 + "TT" = 0, 



(74) 



with H = H/Hq, U = U /Hq and the prime denoting derivative with respect to z, where 

PR = ^, (75) 



PR 



P°R 



3^(1 + z)\ 



H 2 (l + z 



,2^/2 



+ U (Ae a 
- 2U (Ae a(t> - Be~ a 



Be 



-atj> 



(76) 
(77) 



From (1771) . putting the condition p$ = p^ (see Fig. 3) at z = 0, we arrive 
at the condition (f)(0) = and the relation A = B. At z = we also have that 
A/>(0) = 3f2° = 3(1 — 0, R ), which implies 0'(O) = ^6(1 — The Friedmann equation 
for z = furnishes H(0) 2 = Q R + Q^ = 1. These are the initial conditions for the system 
( J73|) -( J74|) . To illustrate a solution, we adopt U = 15, A = B = 0.9 and consider at z = 
the boson and radiation fields as the relevant fields in the Universe, with Q° R = = 0.5. 
Note that z = is not the present time, but the time when the important fields are the 
boson and the radiation. 

In Figure 5 the density parameters of the boson and radiation fields are plotted. 
From this figure we can observe that with an equation of state of the form p^ = 
(kinetic dominant) for the boson field in the presence of the radiation field, if we go 
back in time the radiation can dominate the field 0. Then the present model produces 
the conditions that could favor a big bang nucleosynthesis. 
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7. Conclusions 



Starting from an action that describes a boson, a fermion and a common matter 
field, without specifying their potentials, the Noether condition is applied to the 
corresponding Lagrangian and the undefined potentials are constrained. With the 
Lagrangian satisfying the Noether symmetry, the fermion field presents a term of mass 
and the boson field a potential that can drive an accelerated expansion. 

The first result is that the fermion field has a pressureless matter behavior and 
can play the role of dark matter. Additionally, the existence of a Noether symmetry 
proved to be useful and the complete integration of the field equations was possible. 
The obtained solution is then constrained by observational and physical requirements. 
As a result, the boson field presents an equation of state with its value varying quasi- 
periodically in the range [—1, 1], which means that it behaves as matter as well as dark 
energy. In the beginning of the matter dominated era the boson field plays an important 
role, when it behaves as a matter field. A general consequence of this solution is that 
the model can describe the present accelerated period and predicts a Universe with 
oscillating expansion. Further, the Universe evolves in such a way that in the distant 
future its expansion tends to a constant rate. This is, an eternal accelerated expansion 
as usual does not occur in this model. 
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